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Abstract 
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outgoing gravitons. We find a discrepancy between the amplitude calculated from the 
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1 Introduction 



One of the exciting developments to emerge from the work in recent years on string duality is the 
conjecture that in the infinite momentum frame, string theory (M theory) is described by a large 
N matrix model . In addition, the finite N matrix model has been conjectured to describe 
M-theory on a compact light-like circle [fjj . 

In eleven dimensions, this matrix model is just a supersymmetric quantum mechanics with 16 
supercharges. In dimensions 8 — 10, it is a field theory [^, |], ||. In lower dimensions, the story is 
more complicated [7], ^ . 

The matrix model conjecture has been well tested in processes involving scattering of two 
gravitons |§ and other two-body interactions 0, [II], [L2||. Matrix models also reproduce known 
string dualities in remarkable ways, and seem able to reproduce the light-cone three string vertex 



[13|. Seiberg has recently explained how the conjecture might be derived 



In this note, we describe a further test. We will compute, in this note, an amplitude in 
supergravity involving six gravitons - three incoming, three outgoing - and compare this to a 
calculation in M(atrix) theory. We will work in the limit of small velocities and very small 
momentum transfer, with zero longitudinal momentum exchange. This is analoguous to the well- 
known computation of the v A /r 7 force law between two gravitons as a one-loop computation in 
the matrix model |], [15|, |16j. A rough estimate suggests that the matrix model is on the right 



track. If all three gravitons are separated by a similar distance, R, then the 3^3 amplitude 
should behave as -^4, and this is indeed the behavior one obtains by simple power counting on 
the matrix model diagrams. 

However, when the amplitude is examined in more detail, we seem to find a difficulty We 
will take a limit where one graviton is far from the other two, i.e. r 13 ~ r 23 = R ^> r 12 = r. We 
will expand the amplitude in powers of 1/R. In momentum space this corresponds to a limit in 
which one of the momentum transfers, say qx, is much smaller than the other two (q 2 , q 3 ). 

Since the propagator for a graviton goes as 1/R 7 (there is no longitudinal momentum transfer), 
we expect the leading behaviour to go as -ghpr- The vertices are bilinear in momenta, so we might 
expect a term in the amplitude of the form 

(k x ■ k 2 )(k 2 ■ k 3 )(k 3 ■ kx) 



(i-i) 



where kx,k 2 , k 3 are the momenta of the incoming gravitons. 



In the next section we argue by simple power counting on the appropriate Feynman diagrams 
that such a term can not occur in the matrix model. If the matrix model is truly to reproduce 
the supergravity amplitude, then the coefficient of the term in the supergravity amplitude must 
be zero. 
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The third and fourth sections are devoted to calculating the coefficient of the term in super- 
gravity . We find that the coefficient is nonzero. We first do the calculation using the Feynman 
diagrams in supergravity; as a check, we then perform the calculation in string theory. The two 
computations agree. 

Finally, we speculate on possible origins of the discrepancy we have found. 



2 The matrix argument 



In the matrix model, it is easiest to study scattering with zero p + exchange, and at large 
impact parameter, corresponding to large expectation values for the x's. In momentum space, 
this means one studies momentum transfers small compared to the momenta themselves. If the 
incoming and outgoing momenta of the first graviton are k\ and k 2 , respectively, then k\ = k 2 + qi, 
where |<fi| <C \k±\, q + = 0. Similar remarks apply to the momenta k 3 ,k± and k 5 ,k 6 associated 
with the other lines. 

Among the invariants relevant to this problem are the invariant energies associated with the 
various subsystems, k\ • k 3 , k\ • k$, k 3 • k$. In light cone variables, these have a simple form. For 
example, 

k 1 -k 3 = ktk+(v 1 -v 3 ) 2 . (2.1) 
The terms we are looking for in the amplitude, then have the form 

(fci-fc 2 )(fc 2 -fc 3 )(fc3-fci) = ^^+2^+2) (Vl - V 2 ) 2 (V1 - V3) 2 (V3 ~ ^ ^ 

We want to ask how such a term can arise in the matrix model. Without loss of generality, 
take x 3 = 0, \x 2 \ = r, \x±\ = R. In order to generate a term in the matrix model of the form 
above, we must study loop graphs with the following properties: 

• They must depend on two masses ( R and r) so they must contain at least two loops. 

• They must have six external leg insertions; two for (v i — v 2 ), two for (v i — t^), and two for 
(v 3 - v 2 ). 

• The four factors of v i all attach to "heavy" modes, with mass (frequency) of order R. 
An example of such a graph is shown in Fig. 1. 
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Figure 1: One of the Feynman diagrams of the matrix model. 



Simple power counting shows that this graph (and all of the others with these properties) are 
suppressed by at least nine powers of R. The three propagators each give 1/R 2 ] extracting four 
velocities gives 1/R 4 , and the loop integration gives one power of R. Thus the graph goes as 
1/R 9 . Supersymmetry cancellations might give further suppression, but already this is smaller 
than 1/R 7 . An identical suppression can be found in the background field method used by |I| [J. 



One can alternatively do the counting in terms of effective operators. Integrating out first 
the more massive fields, i.e. integrating over the R loop, yields an effective local operator built 
of the light fields. Since X\ doesn't couple directly to light fields, this operator must be of the 
form {v\ — V2) 2 {v\ — v s ) 2 X 2 . Here X represents operators which can couple to the light fields. On 
dimensional grounds, this should go as 1/R 9 . The integral over the light field loop cannot induce 
compensating powers of R. 

We can also ask what sorts of terms are generated by the iteration of the one loop matrix 
model Hamiltonian. It is easy to show that there are no terms of the correct form here as well. 
We will return to this point in our concluding remarks. One can also easily check that higher 
loop contributions are further suppressed. 

In short, we see that no term of the form 

(fci ■ fc 2 )(fc 2 ■ fc 3 )(fc 3 ■ fci) pi -v 2 ) 2 {vi -V3) 2 (v3-V2) 2 

R 7 r 7 K R 7 r 7 W 

can be generated in the matrix model at any loop order. 



3 Supergravity 

We now turn to the calculation in supergravity. In eleven dimensions, this just means that we 
need to consider graviton exchange between the external gravitons, so the required vertices can 
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Figure 2: The supergravity graphs 



be read off the Einstein lagrangian. The relevant diagrams are shown in Figure 2. 



To evaluate the graphs, we need to know the structure of the supergravity vertices. These are 
conveniently collected in |fl7f| . The three graviton vertex appears in eqn. 2.6 of that paper, the 
four graviton vertex in 2.8. Here we give a brief summary of the computation. 

Work in momentum space. Immediately take the limit g 2 ~ —q 3 ^> qi, where qi are the 
momentum transfers. The leading terms in the matrix model don't involve the fermionic variables 
and thus don't change the polarizations; helicity-changing terms are suppressed by powers of the 
impact parameter (or the momentum transfer, in momentum space). So we look for terms in the 
graviton amplitude 

fci ■ k 3 h ■ k 5 k 3 ■ fc 5 £i ■ g 2 & • £465 • £e ,\ 

<£<& 



The vertices where an incoming graviton emits a single, slightly off shell graviton are partic- 
ularly simple. Calling the incoming and outgoing momenta k\ and k 2 , respectively, and q the 
momentum transfer, and taking the polarization indices to be a/i, [3v, aj, the only relevant terms 
(not suppressed by powers of q) are 

The Pj's denote symmetrization in the three gravitons; the subscript denotes the number of 
distinct permutations. Note that in the limit of interest, the g's are small compared to the 
transverse momenta, so we need to keep in every instance the smallest possible number of g's. In 
the first term, then, the only relevant symmetrization is 1 — ^ 2, but this gives nothing new; in the 
second, this does. In this limit, k\ m — /c 2 (all momenta are defined as flowing into the vertex) so 
we have simply 

V a w hrk^rj^rjap. (3.2) 



Next consider the vertex involving three off shell gravitons. For the kinematic structure of 



equation (3.1), We are only interested in the term involving q 2 
term: 

-2P 3 (gi • q 3 'n a i>'nptrihn) 



g 3 pa — gg. There is only one such 
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where the gravitons 1,2,3 carry indices /j,a, u(3, aj, respectively. The symmetrization produces 
no additional factors. 

Thus for the diagram with the three graviton vertex we obtain 

(m 4 )(£ • 3 (h ■ k 3 )(h ■ k 5 )(k 5 ■ k 3 ) (-^) (3.3) 

The second diagram is the only other relevant one for this structure. We need the four graviton 
vertex appearing here. For such a vertex, label the incoming lines k±, /i, a; k 2 v, f3; 93, a, 7; 54, p, A. 
Because this diagram already has -^r, we can ignore all terms of order q at the vertex. In 
addition, we must again insure that the helicity on the incoming and outgoing graviton lines is 
unchanged. A careful examination of the expression for the four graviton vertex indicates that 
there are only two relevant terms, 

-Pi 2 (ki < jk 2p ri J xVt J ,vVa,p) - ZPnih^h^rj^ri^riax). 

The second term does not obviously preserve helicities. However, if one exchanges 2 — > 4, then 
one obtains 

-Pi2(k la k 2p r]^xV^Va,p) - ^PMh^ki^Vx-yV^Vap)- 

In the first term, exchanging 1 — > 2 gives an inequivalent result (but exchanging 3—^4 does not). 
In the second term, exchanging 1 and 2 similarly gives an inequivalent result. So at the vertex 
we have, using k\ = —k 2 , 

-Zk^pki^-qx^-q^riap. 

One can now combine this with the three graviton vertices to obtain for the diagrams (note that 
in this limit, the four-graviton vertex can attach to either of two of the external lines, giving the 
same result) 

(-2i« 4 )(£ • 3 (*i • h)(h ■ h)(k 5 ■ h) f ( 3 - 4 ) 

\Q1Q3J 

The two contributions do not cancel. 

We have done many further checks on these amplitudes. In the next section, we will show 
that this amplitude agrees with a string computation. It is easy to compute many of the other 
terms in the vertex and to compare these with a string computation as well, and we have done 
several additional tests of this sort. 
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4 Strings 



We now turn to the string calculation. We will work in bosonic string theory. The M-point tensor 
amplitude is given by |L8 



.4 



Mifi-^M^MSl ■■S>M 



1XK 



M-2 



J d 2 zi..d 2 z m — 



, [ (zi - Zj) kii exp 



n (zi - Zj) hij exp 

i>j 



E 
E 



d-d 



2 ( Zi z-j 



ki-Cj 



i>3 ~ Z i)* i>3 Z 3 



E 

i>j 

E 



where 



A 



d 2 Z a d 2 Z b d 2 Z c 



\Z a ~ Z\\ 2 \Z a — Z C \ 2 \Z C — ZA 2 



(4.1) 



(4.2) 



means that three variables are to be given arbitrary values and not integrated. Here we have 
introduced 

kij = ki ■ V 4 (4-3) 
We are interested in terms where the polarizations dot into themselves, thus we need 

C1.C2C1.C2C3.C4C3.C4C5.CeC5.C6 (4.4) 
We will in addition choose z 6 = 00, z 5 = 1, z 4 = 0. This means that all terms involving z 6 cancel. 

The resulting expression has the structure (we will omit the integration symbol henceforth) 

fci2-2 fcl3 fcl4 fcl5 /C23~fc24~fc25~fc34-2 fc3 5 (A r\ 

^12 ^13 *14 *15 ^23 z 24 z 25 z 34 ^35 \^-°) 

multiplied by its complex conjugate . 

We are interested in poles involving To isolate singularities of this sort, the 

following change of variables is helpful: 

zi = (n + l)z 2 z 2 = z 2 z 3 = r 3 z 2 

The integrand then becomes 



(4.6) 



-2+fc56„-2+fci 2T ,-2+fc34 f( r „ \ 

r l r 3 J[ri,r 3 ,Z 2 ) 



(4.7) 



where 



/(n, r 3 , z 2 ) = (1 + n - r 3 ) fel3 (l + n) k ^(l -z 2 - r lZ2 ) k ^(l - r 3 ) fc23 (l - z 2 ) k ^(l - r 3 z 2 ) k * 5 (4.8) 
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In this form, the poles of interest come from the integration regions where the variables go to 
zero. 

Indeed, it is now a simple matter to expand the function / in, say, ri,r%. Keeping the first 
order terms isolates the pole (kwku) 1 - The remaining z 2 integral is readily done exactly. The 
most singular term, proportional to (k^k^k^)" 1 agrees with the supergravity calculation. The 
term which interests us, proportional to k^k^k^kuk^) -1 is non-vanishing and also agrees with 
the supergravity result. A more detailed study of this integral shows that the properly normalized 
amplitude reproduces our supergravity result. Those attempting to verify these statements may 
find the integral: 



helpful. As noted earlier, we have also checked many of the other terms in the amplitude against 
the supergravity computation. 



5 Implications 

We have not found a satisfactory explanation of the discrepancy we have found. However, there 
are several potential issues in the calculation. 

We must emphasize that we have done a matrix model calculation at a finite value of N. 
This means that we have to compare it to M-theory with a compact lightlike circle as proposed 
in . One might worry that this brings in complications while evaluating the supergravity 
amplitude. However, one expects that for a large value of the light-like radius R (or, equivalently, 
a large value of N), the theory on a compact lightlike circle tends to classical 11-dimensional 
supergravity with corrections that are suppressed as j^. Hence we do not expect corrections from 
the lightlike compactifications to alter our supergravity result (as the corrections scale differently 
with R and N.) It is possible that this is too naive and that DLCQ M-theory contains terms 
which are not present in the classical supergravity lagrangian. However, this seems to lead to the 
counterintuitive conclusion that DLCQ M-theory does not behave as classical supergravity in the 
low energy limit. 

It is also possible that we have blundered in some way in our evaluation of the classical 
amplitude. This would be the simplest resolution of the puzzle. But we have performed numerous 
checks of the Feynman rules, as well as of our approach to treating the string computation. The 
detailed agreement of the two computations is impressive. 

Another concern is that the iteration of the one loop matrix model Hamiltonian will generate 




TT(-l) 



Y(A + 1 - ni)r(B + 1 - n 2 )Y(n 3 + n 4 - A - B - l)Y{B + 1 - n 4 ) 
Y{-B)Y{\ + B)T{A + B- ni -n 2 + 2)Y{n 3 - A) 



(4.9) 
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Figure 3: The gravity graph corresponding to the one- loop iteration. 



contributions to the three graviton scattering amplitude. However, these have the wrong N and R 
dependence to resolve the discrepancy. Indeed, the most singular parts of these graphs are easily 
seen to reproduce the Feynman diagrams in which a graviton is first exchanged, say, between the 
first and second line, and then another graviton is exchanged between the second and third. 

Let us first review how the N counting goes for the one loop diagram, corresponding to the 
2^2 process. On the matrix model side, the result is proportional to N j^ 2 {vi—v 2 ) A /t. How does 
this compare to the supergravity calculation? The invariant amplitude behaves as n 2 {k\ ■ k2) 2 /t. 
However, this is for relativistically normalized states. To go to non-relativistic normalization, we 
need to divide by a factor of \fk+ for each external line. This leaves us with 

i^k+k+iv! - v 2 ) 4 = - v 2 ) 4 - 

This is exactly as above. 

Now let's consider the Feynman diagrams in the 3 — > 3 process. First consider the supergravity 
side. The invariant amplitude is proportional to (ki ■ k^) 2 {kz ■ k^) 2 jk 3 ■ q$, where q$ is the small 
momentum transfer, and we have used the fact that k\ = 0. In terms of light cone variables, this 
is 

{ktkj fjv, - V 2 ) 4 (ktkt) 2 (v2 ~ V3) 4 

h ■ Q3 

The external state normalization factors give (kf k^k^)^ 1 . In terms of iV-dependence, this leaves 

N 1 NlN 3 (v 1 -v 2 ) 4 (v 2 -v 3 ) 4 
R 5 h ■ q 3 

It is easy to see how this is reproduced by the matrix model. The iteration of the lowest order 
Hamiltonian reproduces the velocity factors and gives NiN^N^/ R 6 . The energy denominator is 
(£3 • 9*3) / 'N 2 R^ . So we obtain the amplitude above. 
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One could imagine all sorts of corrections to the matrix model Hamiltonian (at one loop) 
which would give things like ■ ks in the numerator, but all contributions will still have the same 
iV-dependence. This is not the iV-dependence of the contribution we want to cancel. (It is also 
difficult to see how one could obtain the correct f-dependence.) 

Another possible concern is that with so many legs, it is not straightforward to work in light 
cone gauge for generic momenta and polarizations. But by studying particular cases one can 
check that this is not a problem. For example, one can take two of the incoming (and outgoing) 
polarizations to be the same. In this case, it is easy to see that there are no new contributions 
with the same dependence on the momenta as those which we have studied here. 



It may be that there is a conceptual problem with our understanding of the matrix model. 
Perhaps there are additional contributions which should be included which can reproduce the 
missing terms. Such phenomena are familar in light-cone field theory where it is necessary 
to include contact terms arising from integrating out backward-moving particles. In previous 
calculations, such terms have been suppressed by powers of N. However, in the present case, 
they potentially have the correct N dependence to generate the missing term. We have not been 
able to understand in any detail how this might work; if it is the resolution, the question will be 
whether there is some simple rule to generate these extra terms in the framework of the matrix 
model. We should note that the arguments of |L4| strongly suggest that this cannot be the source 
of the difficulty. It may also be that the use of the effective Hamiltonian is problematic. For 
example, the effective action presumably contains acceleration terms, which may complicate the 
Hamiltonian treatment. Finally, there is of course the possibility that the finite N matrix model, 
as presently formulated, does not reproduce DLCQ M-theory. 



One might wonder if the derivation by Dijgraff, Verlinde and Verlinde |]I3| of the three string 
vertex already implies that all graviton scattering amplitudes must be reproduced correctly in 
the matrix model. This is not straightforward. Light cone string theory contains an infinite set 
of additional terms (the so-called contact terms) that are not derived from the string vertex, but 
instead must be put in by hand to restore Lorentz invariance Q. One such term is the famous 
v 4 /r 7 term that we have already referred to. The term of interest in this paper is another such 
term. The question is then whether some property such as analyticity of the amplitude implies 
that the contact terms must be reproduced correctly if the string vertices are reproduced (along 



the lines of |19|.) We do not know the answer to this question. One must also note that higher 



order corrections to the three-string vertex which were undetermined in [13] may be important. 
The most impressive and concise argument to date for the validity of the matrix model is 



that of |L4]]. We do not currently understand how this argument might be reconciled with the 



apparent discrepancy we have found. 



There is also an issue as to whether the discrepancy we have found disappears in the infinite 
momentum limit. This can happen because the size of the zero-brane clusters forming a graviton 

4 We are indebted to Steve Shenker for a crucial discussion of these issues. 
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grows with N. In the limit as N tends to infinity, the clusters are infinitely large, and there is 
no sense in which the separation between the clusters can be assumed to be large. This in turn 
means that we cannot find the effective action by integrating out massive strings as we have done 
here. To do the calculation in this limit, one needs more knowledge of the bound state wave 
function. It is hence conceivable that the discrepancy vanishes in the infinite momentum limit. 



Note added: After the completion of this work, we received a paper by M. Douglas and H. 

which describes a discrepancy between the matrix model and supergravity in the 

It seems likely that this discrepancy 



Ooguri [21 



presence of curved backgrounds, following the analysis of jzO 
is related to the one described in this work. 
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